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We compare the correlation times of the Sweeny and Gliozzi dynamics for two-dimensional Ising 
and three-state Potts models, and the three-dimensional Ising model for the simulations in the 
percolation representation. The results are also compared with Swendsen- Wang and Wolff cluster 
dynamics. It is found that Sweeny and Gliozzi dynamics h ave essentially the same dynamical critical 



behavior. Contrary to Gliozzi's claim (cond-mat/0201281: ), the Gliozzi dynamics has critical slowing 



down comparable to that of other cluster methods. For the two-dimensional Ising model, both 
Sweeny and Gliozzi dynamics give good fits to logarithmic size dependences; for two-dimensional 
three-state Potts model, their dynamical critical exponent z is 0.49 ± 0.01; the three-dimensional 
Ising model has z = 0.37 ± 0.02. 

PACS numbers: 05.50.+q, 05.10.Ln, 75.10.Hk 



Cluster algorithms have played an interesting role in 
statistical physics, both for their importance in con- 
structing efficient computational algorithms and their 
unusual dynamical properties. Very recently, Gliozzi has 
introduced a new cluster algorithm, which he claims to 
be free of critical slowing down |j| . Since an under- 
standing of the dynamics at the critical temperature is 
central to both dynamic universality classes and compu- 
tational efficiency, we have made a comparison of the four 
main cluster methods for Potts models due to Sweeny Q , 
Swendsen and Wang [|, Wolff g, and Gliozzi §. 

In 1983 Sweeny Q simulated the Potts model || di- 
rectly in the Fortuin-Kasteleyn representation |^|, given 
by the probability distribution of percolation bond con- 
figuration r, 



P(T) oc 



1-p 



(1) 



where p = 1 — exp(— J/(kT)j, J is coupling constant 
in the Potts model, k is Boltzmann constant, and T is 
temperature; b is the number of bonds present on the 
(hypercubic) lattice, and iV c is the number of clusters of 
the percolation configuration. 

Sweeny used a heat-bath rate of "flipping" the links 
from occupied to empty or vice versa. Consider a par- 
ticular link. We define the state 1 to mean a presence of 
a bond, for absence of a bond, (3 for unoccupied link 
such that the two nearest neighbor sites are on differ- 
ent clusters, and 7 for unoccupied link that spans the 
same cluster, see Fig. [|. In addition, we define lp to be 
an occupied bond, the removal of which leads to /3, and 
similarly 1 7 is an occupied bond, the removal of which 
leads to 7. Finally, a dot • represents an arbitrary state 
(0 or 1). Assuming a link is picked up at random, the 



FIG. 1: Different types of links: 1 (occupied bond) , (absence 
of a bond), j3 (unoccupied link that bridges two clusters), 
7 (unoccupied link that spans the same cluster), lp (bond, 
removing of which breaks the cluster), 1 7 (bond, removing of 
which does not break the cluster). 



transition rate for Sweeny's dynamics is 
w{- -> 1 7 ) =p, w( 



w(-->W = 



P 



(l-p)q+p 



7) = 1 ~P, 

(1-P)« 



(l-p)q+p 



(2a) 
(2b) 



Note that the transition probabilities do not depend on 
the initial state, as it is generally the case for heat-bath 
algorithms. 

Gliozzi also worked in the Fortuin-Kasteleyn represen- 
tation |(| and proposed a variation on the flip rate jl]] 
that is only slightly different from Sweeny's: 

w(l -> 1) = p, ti;(l -► 0) = 1 - p, (3a) 
10(7 — ► 1) = p, w(7 — > 0) = 1 — p, (3b) 
w{f3^l)=p/q, w(P^0) = l-p/q. (3c) 

One can easily show that both of algorithms satisfy 
detailed balance with respect to the equilibrium distri- 
bution, Eq. (nl). Both Sweeny and Gliozzi's rates are 



applicable to model with real positive q (not just integer 
values), although Gliozzi's rate has the restriction that 
p/q < 1. We also note that Chayes and Machta proposed 
a cluster algorithm for noninteger q that takes 0(L d ) op- 
erations per sweep [0. 

A key implementation issue in both algorithms is the 
determination of whether two sites belong to the same 
cluster. For two dimensions, Sweeny devised a very effi- 
cient method based on the special topology. In three and 
higher dimensions, one has to contend with the short- 
comings that a move takes an amount of CPU time pro- 
portional to the typical size of the clusters. 

To separate the issues of the intrinsic dynamics of the 
transition rates from the efficiency of the numerical im- 
plementation, we present results for the correlation times 
of both Sweeny and Gliozzi dynamics, measured in Monte 
Carlo sweeps regardless of the actual computer time for a 
particular implementation. We also compare them in the 
same way with Swendsen-Wang Q and Wolff pj cluster 
dynamics. 

We consider the integrated correlation time, which is 
directly relevant to the magnitude of statistical errors || . 
Instead of the usual method of computing the autocorre- 
lation function for, say, energy, 
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a totally equivalent way is to compute the variance of M 
consecutive terms of the sum of energies, 



M 



We have for the variance 
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A straightforward variance estimator, i.e., the squared 
mean minus the sample mean squared, is used. The usual 
correlation time r is the limit of M going to infinity, and 
is related to the correlation function by 



t = lim tm = 



(7) 



Kikuchi et al have used a similar method to ex- 
tract exponential correlation times. Other conventions 
for specifying the correlation time are in use. For exam- 
ple, Muller-Krumbhaar and Binder M call r' the correla- 
tion time, where r = l + 2r', and Baillie and Coddington 
use the definition r" = r/2. 

The truncated sum converges to its limit according to 
1/M. Thus we obtain the limit with a plot of tm ver- 
sus 1/.M, as shown in Figure |^. The method appears 
excellent for dynamics with relatively small values of r. 
When r is large, we have to go to much bigger values of 
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FIG. 2: tm versus 1/M for the two-dimensional Ising model 
with Gliozzi dynamics on a 32 x 32 lattice, with 8 x 10 6 Monte 
Carlo sweeps. The straight is an error weighted least-square 
fit. The intercept gives the correlation time r = 6.53 ± 0.01. 




FIG. 3: Correlation times r versus lattice linear dimension 
L on a semi-logarithmic scale for the two-dimensional Ising 
model at T c . The circles are for Sweeny, squares for Gliozzi, 
diamonds for Swendsen-Wang, and triangles for Wolff dynam- 
ics. Straight lines are fits to asymptotic slopes at large L. 
The Swendsen-Wang data are from ref. jl3] and Wolff data 
are from ref. |fi"lf . 



M, where the noise can dominate the signal, rendering 
the method less useful. 

We performed our simulations at the critical temper- 
ature kT c (q)/J — l/ln(^/g+ 1) in two dimensions, and 
at J/(kT c ) — 0.221657 for the three-dimensional Ising 
model For each data point, we have spent about 
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TABLE I: Correlation time r for Sweeny and Gliozzi dynam- 
ics at various system sizes. The number in parentheses is 
standard error at the last digits. 



L 


Sweeny rate 


Gliozzi rate 


2D Ising 


2 


2.6121(3) 


2.9414(4) 


4 


3.3280(3) 


3.8722(7) 


8 


4.006(2) 


4.738(1) 


16 


4.688(5) 


5.607(5) 


32 


5.46(2) 


6.51(3) 


64 


6.21(3) 


7.54(5) 


2D q = 3 Potts 


2 


3.0736(2) 


3.8111(4) 


4 


4.5375(8) 


5.9245(5) 


8 


6.424(2) 


8.604(3) 


16 


9.06(2) 


12.18(1) 


32 


12.56(5) 


17.10(5) 


64 


17.5(3) 


24.0(2) 


3D Ising 


2 


3.0242(3) 


3.3456(4) 


4 


4.254(2) 


4.814(16) 


8 


5.594(8) 


6.375(8) 


16 


7.1(1) 


8.0(2) 


32 


9.0(5) 


11.4(10) 



two to four weeks of CPU time on 1GHz Pentium PCs. 
For the Sweeny and Gliozzi dynamics, a unit of time t for 
a lattice linear dimension L in d dimensions is dL d moves 
with bond selected at random. We note that Gliozzi used 
a sequential sweep of the bonds, which gives correlation 
time that is smaller by a factor about 0.6. The correla- 
tion time data are given in Table ||. 

Fig. ^contains a semilogarithmic plot of r versus L for 
the four algorithms: Sweeny, Gliozzi, Swendsen-Wang, 
and Wolff. It is remarkable that Gliozzi dynamics and to 
some extent Sweeny is perfectly logarithmic in size down 
to 2 x 2 lattice. For both Swendsen-Wang and Wolff, 
some curvature can be seen in this plot, and the data can 
also be fitted with a small exponent of 0.25 jll] as well. 
We consider it still to be an open question for Swendsen- 
Wang and Wolff dynamics whether the correlation times 
depend on lattice sizes as a small power or logarithmically 



In Fig. ||, results for the two-dimensional, three-state 
Potts model are presented. The straight lines in these 
log-log plots show good power-law dependence for all 
algorithms. The data show an interesting possibility 
that Sweeny and Gliozzi belong to a different dynami- 
cal universality class from Swendsen-Wang. Linear least- 
squares fits give the dynamical critical exponent z = 
0.56 ±0.01 for Swendsen-Wang and 0.58 ±0.01 for Wolff, 
with lower values of 0.48±0.01 for Sweeny and 0.49±0.01 
for Gliozzi dynamics. The differences are statistically sig- 
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FIG. 4: Same as Fig. |^, but on a double logarithmic scale for 
the two-dimensional three-state Potts model at T c . The num- 
ber on the line indicates the slope of the line. The Swendsen- 
Wang and Wolff data are from ref. Elf . 
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FIG. 5: Correlation times r versus lattice linear dimension L 
on a double logarithmic scale for the three-dimensional Ising 
model at T c of Sweeny (circles), Gliozzi (squares), Swendsen- 
Wang (diamonds) , and Wolff (triangles) dynamics. The num- 
ber on the straight line indicates the slope of the line. 



nificant. 

Figure |^ is for the three-dimensional Ising model. 
Again we see similar behavior, although the data are less 
accurate comparing to that of two dimensions. Sweeny 
and Gliozzi dynamics have approximately the same dy- 
namical critical exponent of 0.37 ±0.02, while that of the 
Swendsen-Wang is larger at about 0.5. Wolff single clus- 
ter algorithm gives remarkably small correlation times 
and a small exponent of 0.29 15 . The Swendsen-Wang 



4 



data show some curvature and might approach the same 
asymptotic value for large lattices. 

In summary, we have computed the correlation times 
for cluster algorithms. It is clear that Sweeny and Gliozzi 
have the same dynamics. In all cases studied here, 
the Sweeny rate is actually better than Gliozzi's. Both 
rates reduce critical slowing down, but not completely 
eliminate it. Both Sweeny and Gliozzi dynamics seem 



to have somewhat smaller dynamical critical exponents 
than Swendsen-Wang when measured in units of Monte 
Carlo sweeps, as we have done in this paper. 
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